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ABSTRACT 

We use recent observations of high-redshift galaxies to study the evolution of 
galactic disks over the redshift range < z 1. The data are inconsistent 
with models in which disks were already assembled at z = 1 and have evolved 
only in luminosity since that time. Assuming that disk properties change 
with redshift as powers of 1 + z and analysing the observations assuming 
an Einstein-de Sitter universe, we find that for given rotation speed, disk 
scalelength decreases with z as (1 + z)^^, total B-band mass-to-light ratio 
decreases with z as ~ (1 + z)^^, and disk luminosity (again in B) depends only 
weakly on z. These scalings are consistent with current data on the evolution of 
disk galaxy abundance as a function of size and luminosity. Both the scalings 
and the abundance evolution are close to the predictions of hierarchical models 
for galaxy formation. If different cosmogonies are compared, the observed 
evolution in disk-size and disk abundance favours a flat low-fJo universe over 
an Einstein-de Sitter universe. 

Key words: galaxies: formation - galaxies: evolution - galaxies: spirals - 
cosmology: theory - dark matter 



1 INTRODUCTION 

The Hubble Space Telescope (HST) and large ground-based telescopes can now provide 
photometric and kinematic data for normal galaxies out to redshifts of about one. Such 
observations are providing new results on the evolution of disk galaxies, showing that the 
population of these objects changes substantially with cosmic time. For example, disks of a 
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given size appear brighter in the past, or equivalently disks of given luminosity appear smaller 
(e.g. Schade et al 1996; Lilly et al 1997 and references therein). A traditional approach to 
interpreting such data has invoked luminosity evolution, assuming that normal giant galaxies 
were already assembled by redshift unity and have evolved primarily in luminosity since that 
time. The purpose of this paper is to examine this hypothesis in the light of recent data on 
disk sizes, luminosities and kinematics. We find it to be untenable if the data are taken at 
face value; strong evolution in the structure of disks appears to be required. The form of this 
evolution is quite similar to that expected in recent models for disk formation in hierarchical 
cosmologies. 

In §2, we use recent local and high-redshift galaxy samples to fit for the values of the 
exponents in power-law representations of the evolution of disk sizes, luminosities, rotation 
speeds and abundances. In §3, we derive the scaling relations predicted for disk properties 
by hierarchical models and we compare them with the observational relations. §4 presents 
some discussion of these results. 

To avoid confusion, we transform all cosmology-dependent observational quantities to 
an assumed Einstein-de Sitter (EdS) universe with h = 0.5 (where h is the present Hubble 
constant in units of lOOkms"^ Mpc^^). We transform theoretical predictions in this way 
also, even when they refer to a different cosmology. 

2 OBSERVATIONAL FACTS 

To study the photometric and kinematic evolution of galaxies, we obviously need both high- 
redshift and local samples. For the high-redshift data, we use primarily the Canada-France 
Redshift Survey (hereafter CFRS, see Lilly et al 1995; Schade et al 1996) and the KECK 
kinematical study of Vogt et al (1997a,b). Ideally, we would like to have a local galaxy 
sample that has both photometric and kinematic information. However, at present this is 
not yet available. We therefore make use of several samples (Kent 1985; Pierce & Tully 1988, 
1992; Courteau 1996, 1997; de Jong & van der Kruit 1994, de Jong 1996). We use simple 
power-law parametrisations of the evolution of observational relations (size-rotation speed, 
size-luminosity, abundance-luminosity, TuUy-Fisher, etc.) and we determine the evolution 
exponents empirically from the data. 
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2.1 The Rd-Vc Relation 

We first study the relation between disk size and rotation speed. The local sample is from 
Courteau (1996, 1997) whereas the high-redshift data are taken from Vogt et al (1997a,b). 
There are respectively 306 and 16 galaxies in these two samples. Figure la shows the his- 
togram of logQrv = log-Rd/Vc, where Rd (in kpc) is the disk scale-length and Vc (in kms~^) 
is the rotation speed. The Rd distribution at fixed Vc is expected to scale in proportion to 
Vc with a coefficient depending on the redshift (see §3); the marginal distribution of Rd/Vc 
should thus isolate the redshift dependence clearly. The dashed line is for the high-redshift 
data while the solid line is for the local sample. Clearly Rd/Vc is smaller at high-redshift, i.e., 
for given Vc, high-redshift disks are smaller (Vogt et al 1997a,b). According to a Kolmogorov- 
Smirnov (KS) test, the probability that these two samples are drawn from the same parent 
distribution is 1.5 x 10~^. Clearly the assumption that disk sizes have not changed since z = 1 
is inconsistent with the data. This excludes models with pure luminosity evolution. To quan- 
tify the size evolution, we adopt simple parametrisation, Qrv{z) = (5rv(0)(l+z)^", and evolve 
each high-redshift galaxy to the present day. We then compare the evolved distribution of 
Rd/Vc to that of the local sample using the KS test. The resulting probability is shown as 
the solid line in Figure Id. One can see that the most probable value of (corresponding to 
the largest KS probability) is = —0.95, with values outside —1.4 < < —0.54 excluded 
at 95% confidence. 

2.2 The Tully-Fisher Relation 

Vogt et al (1997a,b) have studied the evolution of the Tully-Fisher relation with redshift. 
They obtained kinematic data for high- 2; disks using KECK spectroscopy and compared 
them to the local i?-band Tully-Fisher relation of Pierce & Tully (1988, 1992). We have 
converted their quoted data to correspond to an assumed EdS cosmology with h = 0.5 (see 
§3). Unlike the Rd-Vc relation discussed above, the comparison between the high- and low- 
redshift Tully-Fisher relations suffers from the uncertainties in the relative calibration of the 
various photometric systems. We follow Vogt et al in comparing to the Pierce & Tully sample 
even though it contains only 26 galaxies in the Ursa Major cluster, but we note that the 
Tully-Fisher zero-point of this cluster has a significant systematic uncertainty (c.f. Giovanelli 
et al 1997). We define the quantity Qtf = L/Vc, where L denotes galaxy luminosity. Since 
L is roughly proportional to V^, the distribution of Qtf can be used to measure evolution of 



© 1997 RAS, MNRAS 000, 



4 



~n — I — I — I — I — I — I — I — r~| — I — I — I — I — I — I — I — I — r 

:a) 

I 
I 



I I J-TTU-I I I I I M I I lU-i I I I I I 



2.5 
2 

1.5 
1 

0.5 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r~i — I — r 



b) 



_L 



I I I I I I 1 1 



-2.5 -2 -1.5 -1 
Log R,/V, 





-0.5 0.5 1 1.5 2 2.5 
-3 Log(Vj-M,/2.5 




-3 -2 
Log(R,) + M3/7.5 



Figure 1. Histograms for various physical quantities for local (solid lines) and high-redshift (dashed lines) samples. The area 
under all the histograms is one. (a) logQrv = ^og R^/Vc- The local sample is from Courteau (1996, 1997) and has 306 galaxies. 
The high-redshift sample from Vogt et al (1997a, b) has 16 galaxies, (b) logQtf = —31og{Vc) — Mg /2.5. The local sample is from 
Pierce & TuUy (1988) and has 26 galaxies. The high-redshift sample is the same as in Fig. la. (c) log Qrm = log(-Rti) + Mb /7.5. 
The local r-band sample is from Kent (1985) and has 51 galaxies. We have converted it to B using B — r = 0.9 (Kent 1984). 
The high-redshift sample from the CFRS survey (Schade et al 1996) has 103 galaxies. The solid and dashed curves are normal 
fits to the local and high-redshift data, (d) KS test probability as a function of evolution exponents for the three physical 
quantities plotted in panels a, b and c (see text). 

the zero-point of the Tully-Fisher relation. Histograms of logQtf = — Mb /2.5 — 31ogVc are 
shown in Fig. lb. Here Mb is the B-band absolute magnitude of the galaxies. It is clear that 
the data are consistent with no zero-point evolution, in agreement with the conclusion of 
Vogt et al (1997a,b) and Hudson et al (1997). A KS test shows the no-evolution assumption 
to be acceptable. To put limits on evolution of the Tully-Fisher relation, we set Qtf{z) = 
Qtf (0)(l+2;)^*f and perform KS tests on the "evolved" distributions as before. The constraints 
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on /3tf are indicated by the thin dashed hne in Figure Id. With the current data, the evolution 
index is not well constrained, only values outside —1.44 < Pti < 0.36 can be excluded 
with better than 95% confidence. More data, particularly at high-redshift, are needed to 
tie down the evolution of TuUy-Fisher relation, and it is important to ensure consistent 
photometric systems at high- and low-redshift. If we combine the Rd-Vc and Tully-Fisher 
relations, we can estimate the evolution of the total S-band mass-to-light ratio of the visible 
regions of galaxies: M/L oc V^Rd/L = Qr^/Qa oc {1 + 2)^'-v-ftf_ = g, this gives 

M/L oc (1 + z)~^'^. However, if the most likely value of jStf is used, then M/L oc (1 + z)~^'^. 
The exponent here is quite uncertain as a result of the uncertain evolution of the Tully-Fisher 
relation. 

2.3 The Rd-Me Relation 

Next we examine the evolution of Rd-Ms relation. We compare the local (volume-limited) 
sample of Kent (1985) with the high-redshift data of Schade et al (1996). We have checked 
that the independent local (diameter- limited) sample of de Jong (1996) roughly agrees with 
that of Kent, when the different selection functions are properly taken into account. Figure Ic 
shows histograms of the quantity logQrm = log(-Rd) + Mb/7.5 for the local and high-redshift 
samples. Notice that for a tight Tully-Fisher relation of the form L oc V^, Qrm should be 
distributed similarly at each redshift to the Q,.^ discussed in §2.1. For a no-evolution model, 
the high- and low-redshift histograms would coincide. Instead, one clearly sees that disks 
of given luminosity are smaller at earlier times, or equivalently, galaxies of given size are 
brighter. A KS test reveals that the probability that these two samples are drawn from 
the same distribution is less than 4.9 x 10^^. Parametrising evolution of the size-luminosity 
relation as Qrmiz) = Qrm(0)(l +2;)^™, we estimate /?rm as before using a KS test. The result 
is shown by the thick dashed line in Figure Id. The most probable value is /3rm = —1.0 
and values outside —1.35 < /3rm < —0.53 can be excluded with 95% significance. This result 
agrees with Schade et al (1996) who found that the luminosity of disks of given size increases 
approximately as {1 + z)^ in the CFRS. If we use only the big disks as defined in Lilly et 
al (1997), then we obtain /5rm ~ —0.5. This indicates that there might be some differential 
evolution in the Rd — Mb relation for small and big disks. Unfortunately, current data do 
not allow a more detailed analysis of this evolution. 
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2.4 Disk Size Functions 

Schade et al (1996) studied the redshift dependence of the abundance of disk galaxies as a 
function of size. By analogy with the luminosity function this quantity is referred to as the 
"size function" . They found that the abundance of big, bright disks (scalelengths larger than 
5kpc and absolute magnitudes brighter than —20) is similar at 2; ~ 0.3 and at 0.75, whereas 
the abundance of small, bright disks (scalelengths less than 4kpc but Mb < —20) is 5 — 10 
times higher at the higher redshift. A comparison of more recent HST imaging of CFRS 
and LDSS galaxies with the local sample of De Jong et al (1996) confirms and strengthens 
this result (Lilly et al 1997). In addition, if all galaxies to the CFRS limit are included, 
the size functions for low- and high-redshift samples are similar. These authors interpret 
their results as suggesting that the abundances and sizes of big disks have not changed 
significantly since z = 1; the observed evolution reflects changes in disk luminosity caused 
by changes in star formation rate. This interpretation conflicts with our earlier conclusions, 
and in this subsection we test whether the observed size functions are also consistent with 
evolution of the kind inferred above. 

We start with the observed luminosity functions for disk galaxies at various redshifts. 
Using an empirical Gaussian fit to the distributions of log{Rii) + Mb/7.5 in Fig. Ic, and 
assuming that for given Mb disk size varies as Ra oc {1 + z)~^ (see §2.3), we can transform 
the luminosity function into a size function at each redshift. The solid curves in Figure 2 
show the size functions derived from the luminosity functions given by Lilly et al (1995) for 
blue galaxies (which are mostly late-type galaxies) in the CFRS. Results are shown for two 
absolute magnitude limits Mb = —20 and Mb = —18.5: the former is approximately the 
magnitude limit for the high-redshift sample (with median redshift 0.75) whereas the latter 
is that for the low-redshift sample (with median redshift 0.3). We see that the abundance of 
disks with Rd > 4kpc hardly changes with redshift while the abundance of smaller galaxies 
changes substantially. Not surprisingly, this agrees well with the conclusions of Schade et al 
(1996) and Lilly et al (1997) since both our evolution rate and our distribution of Rd for 
given Mb were derived from the CFRS data. 

We now compare these size functions with a model based on the kind of evolution inferred 
in previous sections. Let us approximate the (comoving) abundance of galaxies with rotation 
speed Vc at redshift z as 




(1) 
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Figure 2. The size functions for galactic disks. The solid curves show results derived directly from the CFRS, whereas the 
dashed curves show the model of equation (jsj) for /3rv = —1.0, /3n = 3.0, 7 = —4 and Nt, = 1.0 X 10~'^Mpc~^. The lower pair 
of curves is for galaxies with < —20 a.t z = 0.3; the middle pair for galaxies with < —20 at 2: = 0.75; the upper pair 
for galaxies with Mb < —18.5 at z = 0.3. 

where is some fiducial rotation speed, A^^ is the abundance of galaxies, /5n parameterises 
the abundance evolution, and 7 parameterises the shape of the rotation speed function. Note 
that this is a minimal parametrisation, which can be valid only over limited ranges of Vc 
and z. In practice we are interested in values of Vc similar to those of galaxies. If, as 
observations suggest, disks show a reasonably tight Tully-Fisher relation at each redshift, 
we can use the parametrisation of §2.2 to convert this to a luminosity function, 

. / r \(7-2)/3 

N(L, z)dL = -NJl + zf- oc (1 + z)^"-(^+i)^«/3^(^-2)/3dL, (2) 

3 \LJz) LJz) 
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where L^{z) = (1 + q and L^q is the luminosity corresponding to V* at 2 = 0. 

In §2.1 we found that the mean size- rotation speed relation can be well represented by 
Rd oc Vc{l + z)^" . The scatter in Ini?^ about this relation is approximately normal with an 
rms of (T ~ 0.5 (cf. Fig. Ic). This allows us to calculate the size function for galaxies with 
absolute ^-magnitude brighter than Mum and so rotation speed greater than {Ln^/ L^Y^^Vi,. 
We find 



N{Rd,z)dRd = N^{l + z[ 



,/3n 



Rd 



R.(z) 



lax ^ 

X / — exp 

3 v27r(j 



dRd 

R.{z) 



- (1 + 7)3; 



dx\ (3) 



W = [ln(10)/7.5] (Mn„, - M,) + ln[/?rf//?,(z)], (4) 

where is the absolute magnitude corresponding to L^, and Ri,{z) = -Ro(l + z)^''' is 
the median value of disk scalelength for galaxies with = at redshift z. From the 
solid histogram shown in Fig. la we have (i?o/kpc) ^ 0.03(\4/kms~^). For ^ 
and Rd ^ R-k{,z), the integration in equation (0) depends only weakly on i?^, and the size 
function scales with Rd and z as 

N{Rd, z)dRd oc (1 + (5) 

Thus, the observational result that the comoving number density of large and bright disks 
does not change significantly with redshift requires (3^ ~ PniX + 7)- 

The interpretation of Schade et al (1996) and Lilly et al (1997) corresponds to /3n = 
/3rv = and /5tf = 3 - no evolution in abundance or size but strong evolution in luminosity. 
In contrast, our earlier analysis suggested ~ —1, implying /3n ~ — (I + 7). At luminosities 
around the slope of observed luminosity functions is near —2 which, according to equation 
(I), corresponds to 7 ~ —4. Consistency then requires /?n ~ 3. In Fig. 2 we show the 
predictions of equation (||) for = -1-0, {3^ = 3.0, 7 = -4, = 1.0 x lO^^Mpc"^ and 
A^iim = —20. These agree well both with the size functions derived from the luminosity 
functions of Lilly et al (1995) and with the direct observational estimates of Schade et al 
(1996). The figure also shows predictions for disks at z = 0.3 with Mb < —18.5. The size 
function of this low-redshift sample is similar to that of the high-redshift sample limited 
at Mb = —20 as found by Lilly et al (1997) in their observational samples. We conclude 
that the observed size functions are indeed compatible with a model in which disks of given 
rotation speed have a size which varies as (l + z)^^ and an abundance which varies as (1 + z)^. 
If the redshift dependence of the Tully-Fisher relation is weak, then at fixed luminosity the 
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comoving density in this model evolves approximately as {1 + z)^ (see equation 0). This 
agrees with the behaviour which Lilly et al (1996) inferred from an analysis of the comoving 
i?-band luminosity density. 

3 MODEL SCALINGS 

The analysis of the last section shows that current data are inconsistent with the hypothesis 
that disk galaxies were already assembled at redshift unity and have only evolved in lumi- 
nosity since that time. The current section compares with the predictions of hierarchical 
models for galaxy formation, based on the work reported in Mo, Mao & White (1997). To 
clarify the relevant scalings, it is sufficient to treat galaxies as exponential disks embedded 
in isothermal halos; for simplicity we neglect both the disk mass and the halo core radius. 
The luminosity and size of the disk then scale as 

L oc , n K tt^, Rd oc A———, (6) 
Ttot(2:) H{z) HqH{z) 

where Ttot is the total mass-to-light ratio within the halo virial radius and A is the di- 
mensionless spin parameter. For given circular velocity, disk size depends on redshift only 
through the Hubble constant H{z), while the luminosity has an additional dependence 
through the mass-to-light ratio. If we make the approximations Ttot(-2) = Trtot(0)(l + z)^'^ 
and H{z) = Ho{l + ,2)^'% then from the definitions, Qtt = L/V^ and Qn = Rd/Vc, we have 

QAz) = Qtf(o)(i + ^)-^^-'^N Q„(^) = g,v(o)(i + zy^K (7) 

Thus we see that /5h is equal to /?„ as defined above. In an EdS universe, /?h = 3/2, so this 
simple model predicts that for given circular velocity, the mean size of disks should vary as 
(1 + 2;)"^/^. This is somewhat too steep to be compatible with the data analysed in §2.1 (cf. 
Figure Id). 

In order to compare the data with predictions for other cosmologies, we can either reanal- 
yse the observations assuming the cosmology under consideration, or transform the model 
predictions into those expected when the data are analysed assuming an EdS universe. Here 
we adopt the second approach. Cosmology enters our problem through the Hubble constant, 
the angular size and luminosity distances and the comoving volume. For moderate redshifts 
{0 < z < 1), the ratios of these quantities in a general cosmology to those in an EdS universe 
can be approximated by power-laws in 1 + z: 
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5£l = (1 + iM = ii+ = (1 + (8) 

H{z) dA{z) dv/dz 

where the tilded quantities refer to the EdS cosmology; rf^ is the angular diameter distance 
and dv/dz is the differential comoving volume. The indices ah,(yd and can be approxi- 
mated by 

an = -1.5(1 - «<i = 0.78(1 - 1]^'), = 3-1(1 - f^o''), (9) 

for flat cosmologies with fio + ^ = 1 ? and by 

ah = -0.5(1 - ad = 0.35(1 - fi^"), = 1-2(1 - ^^0'°), (10) 

for open cosmologies with i7o < 1 and A = 0. Using the definitions of Qn, Qrm and Qtf 
together with equations (H) and (||) we have 

These allow the predicted evolution exponents for any cosmological model to be compared 
directly with the observed exponents estimated assuming an EdS universe. For example, a 
fiat cosmology with Qq = 0.3 and A = 0.7 has ah = — 1.0,0;^ = 0.43 and a^, = 1.68. Hence 
/?rv = —0.93 and /?tf = —1-36 — (3r- If there is no evolution in the Tully-Fisher zero-point 
(cf. §2.2), then /5t = —1.36, and P^m = —0.93, in reasonable agreement with the exponents 
inferred from the data in §2. 

The values of (3^ and 7 characterizing the number- density evolution of disks (cf. equation 
|jl[|) can be estimated from the Press-Schechter formalism (Press & Schechter 1974). This 
provides simple formulae for the comoving abundance of haloes as a function of redshift and 
of halo mass Mh or circular velocity Vc (e.g. White & Frenk 1991). For currently popular 
cosmogonies, the slope of the abundance curves for Vc values corresponding to the rotation 
speeds of observed disks implies 7 ~ —3.5 —3.8, while the evolution of the abundance 
at fixed Vc in this range gives (3^ ~ 1-3 1.7 for CDM-like models with f2o = 1 and 
/3n ~ 2.9 —>■ 3.3 for a fiat model with Qq = 0.3 and A = 0.7. The larger (3^ in the latter 
model is due to the volume correction to the EdS universe (cf. equation ^). These values 
are quite similar to those inferred empirically in §2, particularly for the flat model with low 
Qq. Thus the observed evolution both in abundance and in systematic parameter relations 
seems compatible with the expectations of hierarchical formation theories. 
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4 DISCUSSION 

We have shown that current data appear inconsistent with models in which disk galaxies 
are already assembled by redshift one and have evolved only in luminosity since that time. 
Assuming that disk properties vary as powers of 1 + z and analysing them under the as- 
sumption of an EdS universe, we find that for given rotation speed, disk scalelength varies 
as ~ (1 + z)~^, total mass-to-light ratio (in the 5-band) as ~ (1 + z)~^, disk luminosity 
(in 5-band) is almost independent of z, and disk abundance varies as ~ (1 + z)^. These 
scalings agree well with the predictions of hierarchical models for flat low-density universes 
but are marginally weaker than expected in an Einstein-de Sitter universe. Our results for 
the evolution of the size, luminosity and abundance of disk galaxies can be compared to 
those of Bouwens, Broadhurst & Silk (1997) who modeled the counts of faint galaxies as a 
function of apparent magnitude, colour and size in the Hubble Deep Field (Williams et al 
1996). Their preferred model corresponds to jS-^ ~ 3.5 and /5rm ~ —1.5, similar to the values 
derived here. 

In the tests described in this paper, our own preferred model is distinguished from the 
simpler hypothesis proposed by Schade et al (1996) and Lilly et al (1997) only because we 
have included kinematic data. Our results therefore depend critically on the high-redshift ob- 
servations of Vogt et al (1997a,b). The apparent evolution of the size-rotation speed relation 
is inconsistent with pure luminosity evolution, which also requires stronger luminosity evolu- 
tion than can be reconciled with the weak redshift-dependence of the TuUy-Fisher relation. 
The Vogt et al sample is quite small and the comparison of its TuUy-Fisher zero-point to 
that of local samples is uncertain. Further kinematic and photometric data for high-redshift 
galaxies are necessary to confirm the conclusions we have reached. Only with the inclusion 
of kinematic information is it easy to separate evolution in the mass and size of galaxies 
from evolution in their stellar mass-to-light ratio. 
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